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Abstract 

Nanoscale inhomogeneities and impurity clustering are often found to drastically affect the mag- 
netic and transport properties in disordered/diluted systems, giving rise to rich and complex phe- 
nomena. However, the physics of these systems still remains to be explored in more details as 
can be seen from the scarce literature available. We present a detailed theoretical analysis of 
the effects of nanoscale inhomogeneities on the magnetic excitation spectrum in diluted systems. 
The calculations are performed on relatively large systems (up to A^=66'^). It is found that even 
low concentrations of inhomogeneities have drastic effects on both the magnon density of states 
and magnon excitations. These effects become even more pronounced in the case of short ranged 
magnetic interactions between the impurities. In contrast to the increase of critical temperatures 
Tc reported in previous studies, the spin-stiffness D is systematically suppressed in the presence 
of nanoscale inhomogeneities. Moreover D is found to strongly depend on the inhomogeneities' 
concentration, their size, as well as the range of the magnetic interactions. We believe that this 
detailed numerical work could initiate future experimental studies to probe this rich physics with 
the most appropiate tool. Inelastic Neutron Scattering (INS). 

PACS numbers: 75.30.Ds, 73.21.-b, 75.50.Pp 
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Disordered magnetic systems, such as transition metal alloys^"-, diluted magnetic semi- 
conductors (DMSs)^"-, diluted magnetic oxides (DMOs)^'^, and manganites^i"— have contin- 
ued to attract a great deal of attention ever since their respective discoveries. The prospects 
of these materials for novel spintronic devices led to a plethora of work from both the ex- 
perimental as well as the theoretical point of view. For a long time the major focus of 
spintronics research was mostly on homogeneously diluted systems (with magnetic impuri- 
ties distributed randomly on the host lattice), which were believed to give rise to the much 
coveted room-temperature ferromagnetism. For that purpose huge efforts to grow systems 
as "clean" as possible are made. However, in contrast to the common belief, this does not 
systematically appear to be the optimal route to room-temperature ferromagnetism. For 
instance, in homogeneously or well-annealed Mn doped IIl-V DMSs it has been impossible, 
so far, to go beyond the critical temperature of 140 K for 5% Mn^^""— . 

These findings led people to look into other possible avenues in the ultimate quest for 
the room-temperature phenomenon. Transmission electron microscopy (TEM) experiments 
revealed the presence of coherent Mn-rich spherical nanocrystals in (Ga,Mn)Asi^'i^, which 
exhibited a Curie temperature of ~360 K. Similar nanoscale clusters were also reportedly 
observed in other diluted materials like (Ge,Mn)^"— and (Zn,Co)0^, exhibiting critical 
temperatures in excess of 300 K. The nanoscale spinodal decomposition into regions with 
high and low concentration of magnetic ions was speculated to be the possible reason for the 
apparently high Curie temperatures in these compounds. Nanoscale inhomogeneities have 
also often been detected in other families of compounds. In manganites, inhomogeneities 
arise due to the interplay between the charge, spin, orbital, and lattice degrees of freedom, 
leading to the coexistence of metallic and insulating phases. For example, scanning tunneling 
spectroscopy (STS) images of Lai^ajCa^^MnOs revealed a clear phase separation just below 
the critical temperature^. It is widely believed that the phase separation is at the origin of 
the well known colossal magneto-resistance (CMR) effect^*^. In high Tq superconductors 
experimental studies often indicate the crucial role of inhomogeneities^"— . For example 



in Ref. 



26 



it was found that regions with weak superconductivity can persist to higher 
temperatures if bordered by regions of strong superconductivity. The possibility of increasing 
the maximum transition temperature by controlled distribution of the dopants has also been 
suggested. 

Despite the existence of numerous experimental studies, the effects of nanoscale inhomo- 
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geneities in disordered and diluted systems still remain largely unexplored on the theoretical 
front. Because of the large supercells required density functional theory (DFT) based cal- 
culations for inhomogeneous systems are considerably difficult. For the same reason the 
essentially exact Monte Carlo studies are unfortunately incapable to deal with these sys- 
tems. Moreover the crucial importance of both thermal and transverse fluctuations calls 
for the inevitable need of an exact treatment of the disorder effects, which implies a real 
space treatment. Recently using a real-space RPA approach, able to handle considerably 
large system sizes, it has been demonstrated that nanoscale inhomogeneities can lead to a 
dramatic increase of the Curie temperatures^l. The question which naturally arises is how 
impurity clustering affects the spin excitation spectrum? Indeed the spin-wave excitations 
provide valuable insight into the underlying spin dynamics of the system. Inelastic neutron 
scattering is a powerful experimental tool in this context, as the dispersion relation and the 
magnon damping can be measured directly and accurately. In contrast to non-dilute sys- 
tems, such as manganites^'^^^ cobaltites^S'^, multiferroics^"^, pnictides^'^^, etc., one can 
find very few detailed studies^^ devoted to spin dynamics in dilute magnetic systems. Note 
that till now no experimental data has been reported for spin excitations in III-V DMSs, 
such as (Ga,Mn)As. 

The aim of the current work is to provide such a detailed theoretical account of the 
spin excitation spectrum of diluted magnetic systems in the presence of nanoscale inhomo- 
geneities. We show that nanoclusters of magnetic impurities can have drastic effects on the 
magnon density of states (DOS), the dynamical spectral function, as well as the spin-stiffness 
in these systems, when compared to the homogeneously diluted case. In our calculations, 
we have assumed for the lattice a simple cubic structure with periodic boundary conditions. 
In order to avoid additional parameters, the total concentration of impurities in the whole 
system is fixed to x=0.07. The inhomogeneities are assumed to be spherical of radii tq 
and the concentration inside each of these nanospheres is denoted by Xm- For simplicity 
we also fix Xj„=0.8 for all cases considered in this study. We define the concentration of 
nanospheres in the system as Xns = Ns/N, where Ns is the total number of sites included in 
all the nanospheres and = L'^ is the total number of sites. A variable P/v = {xin/x)xns, is 
used to represent the percentage of total impurities contained within the nanospheres. For 
a particular disorder configuration the nanospheres are distributed in a random manner on 
the lattice, the only restriction imposed is to avoid any overlapping between them. 
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We start with a Hamiltonian describing Nimp interacting spins randomly distributed on 
a lattice of sites, which is given by the diluted random Heisenberg model 

Hneis = — ^ JijPiPj^i ■ Sj (1) 

where the sum ij runs over all sites and the random variable Pi is 1 if the site is occupied by an 
impurity, otherwise it is 0. The above Hamiltonian (Eq.l) is treated within the self-consistent 
local random phase approximation (SC-LRPA), which is a semi-analytical approach based 
on finite temperature Green's functions. Within this approach, the retarded Green's func- 
tions are defined as Gl^{u) = j^^ Gl^ifye'^^'dt, where Gi^it) = -ie{t){[S+{t), S- (0)]) , ((...) 
denotes the thermal average). This describes the transverse spin fluctuations and the index 
"c" denotes the disorder configuration. After performing the Tyablicov decoupling of the 
higher-order Green's functions which appear in the equation of motion of G^j{uj) one gets, 

oj — u! + le 

a " 

where (Sj) are the local magnetizations which have to be calculated self-consistently. I^&q''^) 
and I^E'q''^) are respectively the right and left eigenvectors, associated with the same eigen- 
value of the effective Hamiltonian H^^^, whose matrix elements are given by 

I 

At finite T, H^ff is non-Hermitian (real non-symmetric) but it is bi-orthogonal^. Hence 
one needs to define the left and right eigenvectors of the effective Hamiltonian. However, 
since in the present case all calculations are performed at T = K, the matrix H^ff is real 
symmetric and the left and right eigenvectors are identical. For more details on the SC-LRPA 
approach one can see for example Refs. [s 3^ 39. The SC-LRPA was previously successfully 



implemented to calculate the magnetic excitation spectra in the nearest-neighbor diluted 
Heisenberg model^ as well as in the case of optimally annealed (Ga,Mn)As^. In the latter 
case, an especially good agreement with experiments^!^ was also found. 

Now it is well known from first principles based studies as well as model calculations 
that the exchange couplings in III-V DMS compounds are relatively short-ranged and expo- 
nentially decaying in nature^i^i^. Thus, without any loss of generality, we have assumed 
the interactions between the magnetic impurities to be of the exponential form Jjj= Joexp(- 
|rj-rJ|/A), where A is the damping parameter. For this damping parameter, we focus here 



on two particular values, A = a and a/2 (a is the lattice spacing), corresponding to rel- 
atively long-ranged and short-ranged couplings respectively. In fact it was only recently 
shown that just varying A within this scale can give rise to rather spectacular effects on 
the Curie temperatures in these inhomogeneous diluted systems^I. Our current goal is to 
provide a detailed analysis of the effects of nanoscale inhomogeneities, in the presence of 
these interactions, on the dynamical magnetic properties of diluted systems. 

We first start with the calculation of the magnon DOS for the homogeneously diluted 
case. The average magnon DOS is given by pavgi^) = jf—'l2i{Pii^))c, where Pi{ou) is the 
local magnon DOS, which reads pf{uj) = — ^^^g^^ lmG^^{uj) . In what follows {...)c denotes 
the average over disorder configurations. Figure [1] shows pavg as a function of the energy u, 
for the two different values of A mentioned above. The average over disorder is performed 
for a few hundred configurations in all the following calculations, pavg is found to exhibit a 
regular Gaussian-like shape for the case of A=a. The broad peak is located at 0A2W with 
a half-width of about 0.36iy {W is the magnon excitation bandwidth). For longer ranged 
couplings, Pavg remains essentially similar to that of A = a. On the other hand, for short- 
ranged couplings (A = a/2), pavg has a more irregular and richer structure. The peak in pavg 
is now located at much lower energy, 0.06H^, and a clear long tail extending toward the high 
energies with multiple shoulders appears. These additional features result from clusters of 
impurities weakly coupled to the rest of the system. These shoulders become even more 
pronounced for shorter ranged interactions. It is interesting to note that a similar kind of 
magnon DOS was obtained in the case of (Ga,Mn)As^. However, in Ref. |4l| the exchange 
couplings used had been directly calculated from the "l^-J" model. It has been found that 
A = a/3 provides a very good fit for these couplings. 

In order to analyze the effects of inhomogeneities, we calculate in addition to pavg, the 
local magnon DOS inside (pj„) and outside (pout) the nanospheres. Unlike pavg, the local 
DOS contains information on the nature of the magnon states. We consider the particular 
case of nanospheres with fixed radii ro=2a, concentration inside Xin=0.8, and four different 
concentrations of nanospheres Xns- The results for A=a are depicted in Fig. O Let us first 
focus on Pavg- Froui Fig. [2]|^a), we immediately notice that a relatively small concentration 
of inhomogeneities (x„s~0.02) causes a significant change in the magnon DOS. Indeed, in 
comparison to the homogeneous case, the excitations spectrum bandwidth is now doubled, 
and pavg has a bimodal structure, with a broader peak at higher energies. With increasing 



6 




FIG. 1: (Color online) Average magnon DOS corresponding to the homogeneous case for (a) A=a, 
and (b) X=a/2. The average concentration is x=0.07. The x axis is in units of u/W, where W is 
the magnon spectrum bandwidth. l^~4Jo and O.SJq, for X=a and a/2 respectively. The system 
size is A^=56'^. 

Xns, we observe a gradual transfer of weight from the low to high energy peak. The low energy 
peak shifts to smaller energies which is consistent with the decrease in the concentration of 
impurities outside the nanospheres. In order to have a better understanding of the features 
seen in pavg, we now analyze pm and pout- We observe that pm remains unchanged in all cases 
and exhibit a very small weight from to 0.7W. Thus the high energy peak seen in pavg can 
clearly be attributed to the nanocluster modes. A careful analysis of a single isolated cluster 
reveals that the first non-zero eigenmodes are located at 0.714^, which explains the very small 
weight in below this value. Note that in Fig. 2(a), the shaded regions correspond to the 
discrete spectrum of an isolated single nanosphere, which is calculated over a few hundred 
configurations (random position of the impurities inside the nanosphere) . The weak variation 
of Pin with respect to Xns, indicates that the disappearance of the discreteness in p^ (as seen 
in the isolated nanosphere spectrum) results mainly from the interactions between the cluster 
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FIG. 2: (Color online) pa denotes the average magnon DOS {pavg), local magnon DOS inside (pin), 
and outside the nanospheres (pout), for four different Xns- Pn is the percentage of total impurities 
inside the nanospheres. The blue dashed curves represent the homogeneous pavg from Figure 1(a). 
The shaded regions in (a) correspond to the calculated eigenmodes of a single isolated nanosphere. 
The parameters are X=a, ro=2a, and Xin=0-8. The x axis is in units of ui/W, where W^AJq is 
homogeneous magnon bandwidth. 



impurities and those outside. The above discussion of pavg and pj„ explains naturally the 
behavior of pout- In the case of more extended couplings, it is expected that (i) pavg loses 
progressively the bimodal nature, (ii) the pseudo-gap in pm at low energies is filled gradually, 
and (iii) the second peak in pj,„ becomes narrower and shifts to higher energies with respect 
to the spectrum of a single isolated cluster. 

We now discuss the effects of short ranged interactions (A=a/2) on the average and 
local magnon DOS. The results are shown in Fig. [31 The magnon bandwidth increases by 
250% with respect to that of the homogeneous system. As seen before, we observe a clear 
transfer of weight in Pavg, from the low to higher energies with increasing Xns- In contrast to 
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FIG. 3: (Color online) 

Pavgi Pirn and Pout for — o/2, for four different x^is- The blue dashed curves 
represent the homogeneous pavg from Figure 1(b). The shaded regions in (a) correspond to the 
calculated eigenmodes of a single isolated nanosphere. The other parameters being same as in 
Fig. [2j The homogeneous magnon bandwidth is H^~0.8Jo. 

the bimodal nature observed for X=a, pavg now exhibits a long wavy tail, extending toward 
higher energies, shows (i) a clear multiple peak structure now, (ii) is independent of 
Xnsi and (iii) a well defined gap of approximately O.SVF is observed. The reasons for the 
appearance of these multiple peaks in are the enhanced discreteness (larger sub-gaps) of 
the eigenmodes of the single isolated nanosphere and the reduced interactions of the cluster 
impurities with those outside. Concerning pout-, besides a shift to lower energies as seen for 
A=a, we now observe that the peak becomes narrower with increase in Xns- (The latter 
feature was absent for the long ranged couplings). The reason for this is with increasing 
Xnsi the concentration of impurities outside decreases and the effective interactions between 
them become weaker. This effect will be even more pronounced for shorter ranged couplings. 
Even though the couplings are comparable, drastic changes between Fig. |2] and Fig. El shows 
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that X=a corresponds to the intermediate range couphngs and A=a/2 definitely to the short 
range regime. Similar kind of dramatic effects have also been reported on the effects of 
nanoscale inhomogeneities on critical temperatures^. 




FIG. 4: (Color online) Average spectral function A(q, a;) as a function of the energy in the (1 
0) direction for different values of Qx, corresponding to the homogeneous case for (a) A=a, and (b) 
X=a/2. The energy axis (x-axis) is in units of uj/W. The system size is A^=66^. (The intensity of 
the peaks have been multiplied by the factors indicated in the parentheses). 

We propose now to focus on the low energy excitation spectrum in these systems. For 
this purpose we evaluate the dynamical spectral function, which provides deeper insight into 
the underlying spin dynamics. This physical quantity can be directly and accurately probed 
by inelastic neutron scattering (INS) experiments. It is defined by A{q^,u) = {A'^{ci,u))c = 
— <^ ^^^^g^^^ ImG"^(q, u;)^ , where {{S^)) = jf—'^iiSi) is the total average magnetization 
over all spin sites, and G''^(q, w) is the Fourier transform of G1j{uj) given in Eq. 2. The 
averaged dynamical spectral function is evaluated from the following 

A{ci,^) = (Y.KW{^-<)^ (4) 

and 

imp 

(5?) 

where Aj = -jj^iy} temperature dependent local parameter. Note that at T = K, all 

A/s = 1. 
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FIG. 5: (Color online) Average spectral function A(q, w) as a function of the energy in the (1 
0) direction, for four different Xns- Pn indicates the percentage of total impurities inside the 
nanospheres. The parameters are A=a, ro=2a, rcj„=0.8, and A^=52^. The x axis is in units oiu/W. 
(The intensity of the peaks have been multiplied by the factors indicated in the parentheses). 



We start with the average spectral function for the homogeneously diluted systems. 



Figs. 4(a) and 4(b) show yl(q, cj), for X=a and a/2 respectively, as a function of the energy 
for different values of the momentum q in the (1 0) direction. First in both cases, as q 
increases the peaks become broader and more asymmetric with a tail extending toward higher 
energies. Well-defined excitations exist only for relatively small values of the momentum, 
beyond q^a ~ 0.247r no well-defined magnons exist. However, for A=a, the well-defined 
excitations persist up to energy values of about 0.25W, whilst for X=a/2 the excitations 
reaches only up to 0.035H^, where W is the A-dependent magnon spectrum bandwidth. We 
remind that 4Jo and O.SJq, for X=a and a/2, respectively. 

We now proceed further and analyze the effects of nanoscale inhomogeneities on the 
dynamical spectral function. The results for A=a are depicted in Fig. [51 As we increase 
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FIG. 6: (Color online) Average spectral function j4(q, (x)) for the case of X=a/2, with the other 
parameters same as in Fig. [5l 

Xns, there is a broadening in the excitations, accompanied with an increase in asymmetry 
and a shift toward the lower energies is observed. These effects are already pronounced even 
for the lowest concentration of nanospheres. For instance, for q^a ~ 0.127r, the magnon 
energies are 0.09W^, 0.075H^, and 0.067H^, for Xns=0, 0.02, and 0.04, respectively. In order 
to analyze the effects of inhomogeneities on the magnon lifetime for a given q, we define the 
ratio -R(q) = 7(q)/a;(q), where 7(q) is the half-width of the excitations. The excitations are 
well-defined in character only when -R(q) < 1. For the aforesaid Qx, the corresponding i?(q)'s 
are 0.2, 0.33, and 0.66, for x„s=0, 0.02, and 0.04, which corresponds to an increase of about 
60% and 200% respectively, compared to the homogeneous case, for these two values of a;^^. 
It is interesting to note that these effects could hardly be anticipated from the magnon DOS 
results (Fig. 2). In fact the analyses of the DOS suggested that the low energy excitations 
should be weakly affected by the inhomogeneities. In the following, we discuss the spectral 
function in the presence of short-ranged interactions, shown in Fig. El As in the previous 
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case, well-defined excitations exist only for small values of the momentum. However, here 
we find that the shift toward the lower energies is strongly enhanced. If we consider the 
particular case of q^a ~ 0.127r, the magnon energies are shifted by 30% and 60% respectively, 
for x„s=0.02 and 0.04, with respect to that of the homogeneous case. The -R(q)'s for this 
value of are 0.4, 0.8, and 1.3 for x„s=0, 0.02, and 0.04. This indicates that the excitations 
have dramatically lost their well defined character as compared to the previous case (Fig. [5]). 
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FIG. 7: (Color online) Magnon dispersion as a function of Qx/t^, for different Xns (0, 0.02, 0.04, 
0.06, and 0.08) for (a) A=a, and (b) X=a/2. The vertical arrows indicate the direction of increasing 
Xns- The solid (red) curve corresponds to the homogeneous case. The shaded regions indicate the 
q values beyond which there is no clear maxima in the magnon excitations, (c) and (d) show the 
normalized spin-stiffness D„ as a function of Xns, corresponding to \=a, and A=a/2, respectively, 
for three different vq: \/^a, 2a, and ^/2a. The concentration inside nanospheres is Xin=0.8. 



From the A(q, uj) calculations, we now extract both the magnon dispersion w(q) and the 
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spin-stiffness We remind tliat in tlie long wavelengtli limit, uj{c\)=D^"'^q^ . In order to 

gauge the effects of the inhomogeneities, we define the normalized spin-stiffness coefficient 
Dn=D^^^ / D^°^ , where _D*"'^ denotes the spin-stiffness of the inhomogeneous systems, and 
jjhom Qf homogeneously diluted system. The values of the spin-stiffness D^°^ are 
found to be 2.9Jo and O.OTJq, for \=a and a/2 respectively, for x=0.07. We would like to 
stress that both u;(q) and D*"'' were obtained from the first non-zero q excitations (q=(f^, 



0, 0)) corresponding to system sizes = L , where L varies from 36 to 52. Figs. 7(a) 



and 7(b) show w(q) as a function of q for various Xns, corresponding to A=a, and a/2 
respectively. First, in both cases, we observe that the inhomogeneities suppress the magnon 
dispersion. This suppression is even more drastic in the case of short-ranged interactions. 
In addition, the quadratic behavior persists up to large q for A=a/2, whilst a clear softening 
of the magnon modes can be seen in the case of the extended couplings, for Qx/t^ >0.2. This 
is consistent with the discussions of Figs. [5] and [61 The well-defined magnon region shrinks 
significantly on going from \=a to A=a/2. This tendency is enhanced on increasing Xns- 
In order to evaluate more quantitatively the effects of the inhomogeneities, the normalized 



spin-stiffness Dn is plotted in Figs. 7(c) and 7(d), as a function of for various radii of 
the nanospheres. As anticipated from the previous figures, one observes for both \=a and 
a/2, a monotonous decrease of the spin-stiffness with increasing Xns- Dn decreases almost 
linearly up to Xns ~0.04. The slope for A=a/2 is twice that of A=a. We find that the 
suppression of the spin-stiffness is greater for larger radii of the nanospheres. For instance, 
for the particular Xns=0.03, the spin-stiffness is reduced by 15% for rQ=\/2a, and by almost 
30% for ro=A/6a, for A=a (Fig. 7(c)), although the concentration outside the nanospheres 
is almost the same. In the case of the short-ranged couplings this effect is even stronger. 
The effects of the nanoscale inhomogeneities on the spin-stiffness are in striking contrast to 
that observed on the critical temperatures, where an increase by more than one order of 
magnitude was found^^. 

To conclude, in this work we have shed light on the effects of nanoscale inhomogeneities 
on the magnetic excitation spectrum of dilute systems. To our knowledge, no such detailed 
numerical study has been performed so far. Compared to homogeneously diluted systems, 
even low concentrations of nanoscale inhomogeneities is sufficient to induce drastic changes 
in the magnon DOS. The magnon dispersion is found to vary significantly with the concen- 
tration of inhomogeneities, leading to a strong suppression of the spin-stiffness. This is in 
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contrast to what has been reported about the critical temperatures. We also find a strong 
increase of the half- width of the magnon excitations in the long wavelength limit. These 
features are strongly enhanced in the case of short-ranged couplings. Now it would be in- 
teresting to study the effects of temperature on the spin dynamics of these inhomogeneous 
systems^. We expect the temperature dependence of the spin-stiffness, especially in systems 
with short-ranged interactions, to be unconventional compared to that of the homogeneous 
systems. Finally, we believe this detailed study would serve to motivate future experimental 
works on these inhomogeneous compounds. 
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